Decoherence, Autler-Townes effect, and dark states in two-tone driving of a 

three-level superconducting system 
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We present a detailed theoretical analysis of a multi-level quantum system coupled to two radiation 
fields and subject to decoherence. We concentrate on an effect known from quantum optics as Autler- 
Townes splitting, which has been recently demonstrated experimentally [M. A. Sillanpaa et at, Phys. 
Rev. Lett. 103, 193601 (2009)] in a superconducting phase qubit. In the three-level approximation, 
we derive analytical solutions and describe how they can be used to extract the decoherence rates 
and to account for the measurement data. Better agreement with the experiment can be obtained 
by extending this model to five levels. Finally, we investigate the stationary states created in the 
experiment and show that their structure is close to that of dark states. 
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I. INTRODUCTION 

During the last decade, the field of superconducting 
qubits has advanced tremendousljji. These systems be- 
have quantum-mechanically, and can be regarded as tun- 
able artificial atoms. Different from their natural coun- 
terparts, they couple strongly with the environment and 
as a result they have shorter coherence times. On the 
other hand, stronger coupling has its own advantages, 
such as easier addressability and shorter gate times. 

A number of quantum phenomena become manifest 
when atoms interact with electromagnetic radiation. 
Among these, electromagnetically induced transparency 
in atomic media allows spectacular effects such as the 
reduction of the group velocity of light to a few meters 
per second^ or even a complete stop. This phenomenon 
requires two fields (the "probe" and the "coupling" - or 
"pump"- field) and a three- level atom. The two fields 
are typically on-resonance (or not far off-resonance) from 
the two transitions. This effect can be used for co- 
herent storage of optical informatior^^, for the realiza- 
tion of quantum repeaters^, for enhanced photon-photon 
interactions^, and for setting up table-top cosmologi- 
cal experiments such as the creation of event horizons^. 
These phenomena are not restricted to atomic physics, 
but solid-state systems can also be used: for example, ul- 
traslow light propagation has been already demonstrated 
in crystals^. Electromagnetically induced transparency 
has its roots in the destructive interference of proba- 
bility amplitudes of the state that is coupled to both 
fields: as a result, the "probe" field is no longer absorbed 
and the atom becomes "trapped" in a "dark" state (a 
superposition formed by the other two states), a phe- 
nomenon called " coherent population trapping" . This 
phenomenon could find applications as a way to precisely 
prepare initial states in future quantum computers. In 
solid-state systems, this effect has so far been demon- 
strated with donor-bound spins in GaAs^, with nitrogen- 
vacancy centers in diamond^, and with single spins in 



quantum dotsiS. 

Superconducting quantum circuits can be operated as 
well as three-level systems interacting with two radiation 
fields on-resonant with the two transitions. So far, this 
has been realized with phase qubits^i'-i2, transmons^^,, 
and flux qubits^'*. From the spectra obtained in these 
experiments it is sometimes difficult to distinguish be- 
tween electromagnetically induced transparency (whose 
origin is quantum interference) and the Autler-Townes 
effeclji^ (which is due to the shift of the resonance of 
the transition frequency which is probed). As a rule of 
thumb to separate the two effectsi^, if the coupling-field 
Rabi frequency (denoted by fJc in this paper) is much 
smaller than the spectral linewidth of the darkened tran- 
sition and in absorption spectroscopy we see a sharp dip 
formed within the linewidth, then we have electromag- 
netically induced transparency. If the Rabi frequency 
of the coupling field is much larger then the linewidth, 
we have the Autler-Townes effect, which appears as the 
splitting of the spectral line into two lines of width sim- 
ilar to the initial one and separated by ~ Oc- Evidence 
for coherent population trapping and electromagnetically 
induced transparency is so far only indirect: there is yet 
no quantum tomography experiment performed to truly 
measure the superposition mentioned above (see the pre- 
vious paragraph). In the experiment of Ref. [TTj an- 
alyzed in this paper the power of the coupling field is 
such that is larger than the linewidth, making the 
second effect dominant. We will therefore refer to the 
phenomenon described below as the Autler-Townes ef- 
fect in a superconducting qubit. 

We consider here a phase qubit: a system consisting 
of a superconducting loop interrupted by a Josephson 
junction^^. The three lowest energy levels are denoted 
by |0), |1), and |2), and the corresponding transitions can 
be driven^^ by microwave radiation fields with frequen- 
cies of the order of a few gigahertz. In atomic physicaiS, 
the states typically form a A configuration, with |0) and 
1 2) being metastable hyperfine or Zeeman levels, while 
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the state |1) is usually an excited electronic state that 
decays at a faster rate. The situation in the case of 
phase qubits differs: since the direct transition |0) |2) 
is suppressed in these systems due to relatively low an- 
harmonicity, the most straightforward operation of the 
phase qubits is as ladder systems, that is, driving only 
on the allowed single-photon transitions |0) — >■ |1) and 
|1) — >■ |2). It is also possible to operate these systems 
in the A configurationi^, by driving two-photon virtual 
transitions from |0) to |2) with an intense microwave field 
at the frequency /02/2. This frequency can be sufficiently 
detuned from the transition /oi that the direct excita- 
tion into the state |1) is small. Depending on the sample 
design, one can measure either the level population of 
the qubit (if it can be addressed directhO, or the scat- 
tered radiation. For example, in Ref. |14l | the qubit is 
embedded in a one-dimensional transmission line, and 
what has been measured is the absorption of a probe mi- 
crowave signal in resonance with the |0) — |1) transition. 
What is seen is that this quantity becomes large (close 
to 1 within a few percentage points) when a more intense 
field is applied to the transition |1) — > |2). 

Future quantum-information processing devices based 
on superconducting circuit architectures could make use 
of microwave-controlled states in three-level systems. In 
addition to the applications listed above, specific to 
the field of circuit quantum electrodynamics one could 
employ three-level effects for the coherent and tunable 
excitation of microwaves in coplanar waveguides--, as 
an alternative way to measure the qubit 's decoherence 
ratesSij for coolingSS and single-photon generations^, as 
a single-atom media quantum amplifier—, and as quan- 
tum switches^^. 

This work is organized as follows: in Section |ll] we 
introduce the phase qubit and derive its corresponding 
Hamiltonian. In Section IIIII we write the general form 
of the Lindblad master equation for the phase qubit ap- 
proximated as a three-level system, including relaxation 
and dephasing in our model. A full derivation of the re- 
laxation term and of the dephasing term can be found in 
Appendices El and [B] respectively. The three- level model 
can be easily extended to account for spurious excitations 
on the next two levels, and the resulting five-level model 
is briefly presented in Section |lVl In Section |V] we then 
explain how to extract the dephasing rates correspond- 
ing to the three levels from the occupation probabilities 
under continuous irradiation (the relaxation rates are ob- 
tained from independent pulsed measurements). Next, 
we present a number of approximate analytical results 
(Section IVip for the Autler-Townes splitting and for the 
occupation probabilities. In Section fVIII we present the 
experimental data for the Autler-Townes effect and we 
compare them with the theoretical predictions based on 
the parameters determined in Section |Vl We show that 
the approximate analytical results presented in Section 
IVII provide a good fit to the data. Better agreement can 
be obtained with numerical simulations for the full five- 
level model of Section ITVl Finally, in Section lVIIII we cal- 



culate the fidelities of the steady states relative to ideal 
dark states. Details about the structure of the steady 
states are delegated to Appendix O 



II. A THREE-LEVEL ARTIFICIAL ATOM 

In this section we introduce our phase qubit^^, and 
analyze it theoretically in the three-level approximation. 
The phase qubit can be pictured as an rf-SQUID (see Fig. 
[IJ, consisting of a single Josephson junction with critical 
current and capacitance C, which has been inserted 
into a superconducting loop with inductance L. We de- 
note by <i>o the fiux quantum, $0 = h/2e = 2.067 x 10~^^ 
Wb. The Josephson energy of the junction is then 
E,j = ($o/27r)Jc; the apphcation of a phase difference ip 
requires an energy Ej{l — coscp), and can be regarded as 
resulting from a nonlinear inductance Lj(ip) = Lj/ cos ip, 
where Lj = ($o/27r)VSj = $o/(27r/e) can be under- 
stood as the Josephson inductance in the limit of small 
phase differences. 
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FIG. 1: Schematic of the phase qubit, consisting of a single 
junction with capacitance C and Josephson inductance Lj in- 
serted into a superconducting loop of inductance L. A nearby 
dc-SQUID is used for single-shot readout and an on-chip ex- 
ternal line is employed to bias the qubit with external flux 
'I'Gxt with both dc and rf components, <l>cxt ~ ^dc + $rf(i). 

The classical Lagrangian of the rf-SQUID can be writ- 
ten in the standard form as kinetic energy minus poten- 
tial energy, 
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($ 



Ej cos 27r — 



(1) 



where <l>/$o is the superconducting phase difference 
across the junction and ^cxt = $dc + ^rf{t) is an ex- 
ternal magnetic flux with both dc and rf components 
applied though an on-chip flux bias coil. By using the 
Legendre transformation, we obtain the Hamiltonian 



2C 2L 

^ ($-$dc)f 
^ 2C 2L 



cos I 27r — 
$0 



(2) 



— Ej cos 2tt — 
$0 
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where the canonically conjugate variable of the fliix $ 
is the charge Q = dC/d^ = C$ accumulated on the 
capacitor C, and we assume \'^dc\ >• |'i'r/(OI- Let us 
first study the time- independent part H' of -ff, defined 
hy H — H' — $$rf(0/^; W6 will return to the time- 
dependent part toward the end of this section. H' can 
be conveniently written as H' = Q^/2C 
the potential energy is 



E'pot, where 



pot 



El 

— ((/'- (/-dc)^ - Ejcosf 



(3) 



with the definitions El = ($o/27r)2/2i, = 27r<I>/<I>o 
and (f>dc = 27r$dc/<&o- Also, the loop inductance L > Lj 
is chosen so that local minima are formed in the potential. 

In order to analyze the quantum states trapped in one 
local minimum of the potential given by Eq. ([3]), we 
expand Eq. (jSj as a Taylor series around the first local 
minimum (pm (implying that < 0m ^ 27r): 



E, 



pot 



Here 



Epoti(p„.) + E'^^,{^„,)A4> + \E';ot{ 



(4) 



Epoti 
E'^oti 

E: 



pot V 



and the first three derivatives are 

= EL{4'ni - 0dc) + Ej sine/)™, 
= El+ Ej cos 0m, 
= -£'jsin0m. 



To obtain the local minimum 0„i , we need to solve the 
equation Epgf{(f>m) — 0. This is not solvable analytically, 
but the quadratic expansion of the sine function around 
37r/2, 



sm( 



1 



(5) 



gives a good-enough approximation of the first positive- 
value local minimum 0™, 




El 



Stt) 



37^)/3-l+r^(6) 



where ^ = (2^/$o)(i//c) = i/i^j = ^^jZ-^l. Since </>™ 
should be a real number, 2 + (20dc - 5n)(3~^ + f3~'^ > 
must be satisfied, thus 



Pdc > 



Stt 



1 

2^ 



/3. 



(7) 



Having found the local minimum 0,„, we can then 
rewrite Eq. Q explicitly, using also cos{(pm) — 
cos{ipm — 3tt/2 + 37r/2) w cos(i^m — 37r/2) cos(37r/2) — 



sin(37r/2) sin(iy9^ 
Spot(A0) « 



37r/2) w (p„ 



Pdc) 



37r/2, 

'Sj(</>m-3V2) 



+ Y ^2 + (20dc - 3^)/3-i + /3-2A02 



+ -^[2^2 + (20dc - 37r)/3-i -f /3- 
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-(20dc-37r + 2/3-i)]A03 



(8) 



The comparison between this cubic potential and the 
original potential is shown in Fig. [5] We see that around 
the local minimum, the cubic potential approximation 
is excellent, which allows us to use it for extracting the 
energy levels. 




FIG. 2: (color online). Potential energy of the rf-SQUID for 
typical parameters Lj — 276 pH, L = 690 pH, and (p^c = 
0.7667r. The minimum is located at (j)m ~ 1.55l7r. The (red) 
solid curve is the exact potential energy calculated with Eq. 
Q; the (blue) dashed curve is the cubic approximation given 
by Eq. (O. The inset shows the energies of the first three 
metastable states |0), and |2) formed in the well centered 
around the local minimum 0™. 

The last step is to derive the quantum Hamiltonian in 
the three-level approximation. To simplify the notations, 
we can define 



A= [2^ 



and we rewrite Eq. 
of A0) 



(20dc 
2(0dc 



37r)/3-i 



13- 



(9) 



as (dropping terms independent 



E, 



pot 



%^A02 



ELi2X + 
12 



A0^ 



(10) 



As one can see from the expression above, the energies 
that define the shape of the metastable well in this ap- 
proximation are EjX and El{2X + £,)■ Also, since the 
value (pm will play no role from now on, we will make 
the convention that the variables (p and = are 
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measured from 4>m as a reference point. With this con- 
vention, we can rewrite Eq. (jlOp in terms of the fluxes 



pot 



1 

2Z* 



1 



(11) 



where L} = Lj/X and L* = 3L/7r(2A + 0- 

Summing up the results, the total time-independent 
Hamiltonian H' can be written as 



where 



2C 



$2 
2L 



(()3 



J 



2$ni 



$3 



2$ni 



Q2 ^2 
~ 2C + 2Z} 



(12) 



(13) 



describes a harmonic oscillator. 

We can now go further and quantize our Hamilto- 
nian. We immediately notice that the canonically conju- 
gate variable of Q is $. Thus the commutation relation 
Q] = ih can be postulated. We can now introduce 
the creation and annihilation operators for the harmonic 
oscillator Eq. (IT51) . and a, defined by 



-(a^+a), Q 



2CuJo 

and satisfying the algebra [a, a] 



hCujQ 



(at 



, (14) 

0, [a,at] = l 

(here luq = y^l/L*jC is the Josephson plasma frequency). 
Thus 



Ho = huo {a^a + 1/2) , 



(15) 



with eigenvectors |n)o, HQ\n)o = hujQ{n + l/2)|n)o. 

To proceed with the full time-independent Hamil- 
tonian H' , we note that around the local minimum, 
$/$o ^ 1- Therefore, we can treat the second term 
on the right-hand side of Eq. ([TT]) and Eq. (fT^ as a 
perturbation. 

Truncating to the lowest three unperturbed eigenstates 
|n)o of Hq we get Eq. (|12l) in the matrix form 



H' = 



hujQ/2 3ry 
377 3fia;o/2 Gy/^j] 
eV^f] 5fiwo/2 



(16) 



where we defined 77 = {h/2Cu!of/'^ /2^oL* . The per- 
thus couples only the states |0)o and 



turbation 

|l)o with strength 3??, and |l)o and |2)o, with strength 
6^21]. The matrix Eq. can be diagonalized ex- 

actly. The resulting eigenvalues are hujQ/2 — llrj'^ /hujQ, 
3hujo/2 - 7hf/huJo, and 5huJo/2 - IQhf/hajQ. For the 
fabrication parameters of our sample (see e.g. Fig. [5]), rj 
is only about 2% of huto; hence the eigenstates \n) of H' 
are quite close to the unperturbed eigenstates |n)o. 

We now turn to the full Hamiltonian Eq. ([2]) , including 
the time-dependent term $$i.f/L, and to avoid carrying 



over the zero-point energy in the equations, we make the 
convention that all the energies are measured from the 
ground-state level. We first give the numerical form of 
the Hamiltonian suitable for the analysis of our experi- 
ments, which is obtained by calculating the matrix ele- 
ments of the time-dependent terms in the basis \n) with 
77 « 0.02/iwo; we obtain 



H{t) = H' 



$$rf(0 

L 

0.695(<) 0.02.g(i) 
0.695(t) c^io g{t) 
0.025(t) g{t) 



where 



m = 



*rt(i) 



(17) 



(18) 



is a coupling constant (measured in hertz), and the tran- 
sition frequencies are 



wio = 0.976ti;o7 ^^21 = 0.952wo 



(19) 



The frequencies wio and ^21 can be found spectroscop- 
ically. For our sample we find wio = 27r x 8.135 GHz 
and UJ21 = 27r X 7.975 GHz. The experimentally de- 
termined relative anharmonicity (wio — W2i)/(a;io + W21) 
is 9.9 X 10"'^, in reasonable agreement with the value 
12 X IQ-^ obtained by using Eqs. (jl9p . 

As a comparison, if we completely neglect the pertur- 
bation (7/ — >■ 0), Eq. (fT7|) simplifies to 



Ho-^^n 



j^git) 

755(i) Wo git) 







g{t) 2u}o 



(20) 



This expression clearly reveals the structure of the circuit 
Hamiltonian as a ladder system, with interlevel coupling 
strengths given by g{t)/^/2 for the |0) — >■ |1) transition 
and g{t) for the |1) — >■ |2) transition. Also from Eq. ((20)) 
as well as from Eq. PT|) one sees that if the qubit is 
weakly anharmonic, the transition |0) — > |2) cannot be 
driven directly. However, at high enough powers, two- 
photon transitions between these two levels are observed 
experimentally. 

In our Autler-Townes splitting experiments, $r/(0 
consists of two microwave tones: the weak probe tone 
$pCos(ajpt) and the strong coupling tone $c cos(ciJci). 
In Eq. (IT71) we can write g{t) as g{t) — gcCos{ujct) + 
gpCos{ujpt), where the coupling and probe amplitudes 
read 



L 



1 



9p = 



1 



L V Chwr, 



The Hamiltonian Eq. (|17p clearly shows that both mi- 
crowave tones are coupled to |0) O |1), |1) O |2) and 
|0) O 1 2) transitions (cross couplings). The cross cou- 
plings introduce additional complications to any analyti- 
cal study of this system. In the next section we will show 
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that they can be ehminated by changing to suitable ro- 
tating reference frames. This aUows us to then derive 
relatively simple analytical results. 



H — ^ (2cti2Ps0'21 — a'22PS ~ PSf22) 
+K{(Joipsa2i + cri2psaio) , (22) 



III. THE MASTER EQUATION 

The driven three-level qubit is inevitably coupled to 
the external electromagnetic environment, which causes 
relaxation and dephasing of the system. One can easily 
extend the well known spin-boson mode P^'^^ of dissipa- 
tion in two-level systems to this three-level case. See 
also Ref. [1^ for a similar treatment of decoherence in a 
multilevel qubit. 

The standard Markovian master equation for the 
Schrodinger-picture density matrix pg of the phase qubit 
has the form 



where k = y^TioT^. Here, (Jij = and the interlevel 
relaxation rates between |1) — ?> |0) and |2) — >■ |1) are de- 
noted as Fio and r2i, respectively. Note also that for a 
three-level system, C-^ci is not just a simple generalization 
of the well-known expression for two-level systems, and 
a mixing term with effective decay constant k appears - 
see the last term in Eq. 1^^. This term however can be 
neglected after performing a rotating wave approxima- 
tion (see below) . A complete derivation of the relaxation 
Liouvillean is given in Appendix El 

For the pure dephasing part of the Liouvillean we will 
use the general expression 



PS 



-j:[H,Ps] 



(21) 



where H is given by Eq. (flT)) . and £ is the total Liou- 
ville superoperator. There are two contributions to C: 
the first is due to relaxation, and we will call it £rci, and 
the second is caused by fluctuations of the energy levels, 
which we will refer to as pure dephasing, £dcp- Thus 
C = £roi + >Cdop- The reason for calling £dcp "pure" is 
that, as we will see below, relaxation produces interlevel 
dephasing as well. A microscopic derivation of these Li- 
ouvilleans is given in Appendix El and [B] One important 
difference between relaxation and pure dephasing is that 
the first results from energy exchange with the environ- 
ment, while the second is caused by energy-conserving 
virtual processes. This allows us to determine the pa- 
rameters entering £rci hi independent experiments, by 
exciting the system and measuring how long it takes to 
decay into the environment (more details in Section IVT) . 

The action of the Liouville superoperator for relaxation 
£101 is given by 



£rcl [ps 



(2(Toi/Oso'io — cTii/Os — Ps^ii) 



£dcp[Ps] = - X! ^^JjPSO-fefc, 



(23) 



where 7^^. — 7^^ are pure interlevel dephasing rates, and 
j,k e {0,1,2}. A derivation of this expression can be 
found in Appendix IB] 

It is convenient to bring the system into a doubly- 
rotating reference frame defined by the operator 



1 

e-'": 









(24) 



Under this unitary transformation the master equation 
preserves its structure, p — —{i/h)[H,p] +£[/?], where 
p — W psU, the Hamiltonian transforms as H — WHU+ 
ih{dW /dt)U, and the Liouvillean as £[p] ~ W£[ps]U. 
Performing a rotating wave approximation (RWA) to 
drop terms oscillating with frequencies cOp, Uc and Wp-I-Wc, 
we obtain 



^(RWA)/^ : 



0.69 (5p + 5ee''^*) /2 



0.69 {gp 



'iSt\ 



/2 



A. 



{gpe-''' + 9c) /2 



(5pe 
A 



iSt 



^9c 
A, 



/2 



(25) 



where the detunings Ap = ujiq — LUp, Ac = W21 — ^c, and for the effective Hamiltonian: 



5 ~ (jjp — uj, 



The oscillating terms in Eq. (j25l) average 
out to zero on the time scale of our experiment, thus 
for the calculation of steady state [p = 0) we can neglect 
them^, thereby obtaining the final matrix representation 



r2p/2 
VLpl2 Ap 17c/2 
Ap + Ac 



(26) 



where the probe and the coupling Rabi frequencies are 
defined as Q.p = 0.69(7p and ilc = 9c, respectively. 
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The relaxation Liouvillean in the doubly- rotating 
frame reads 



i^vc\[p\ = — (2o-oipo-io - CTiip - pail) 

H Y (2o'l2/30-21 - 0-22/3 " PCr22) 

+ K (e"*'^Voipcr21 + e'**CTi2pCTlo) , 

(27) 

and again we will neglect the oscillating term. Also, it 
is easy to check that the form Eq. (^5]) of the pure de- 



Another form for the dephasing Liouvillean used some- 
times in the literatures^ can be obtained from Eq. ([25)1 
by separating the diagonal and off-diagonal terms, 

= r2iP22(o-ii - (T22) + riopii(croo - CTii) 

E ^-fPok^Jk, (30) 

j,fce{0,l,2};jyfe 

where now 710 = 701 = Fio -I- 7^0, 720 = I02 = ^21 + 720 ' 
and 721 = 712 — ^21 + Tio + 721 have the meaning of off- 
diagonal decay rates for the corresponding density matrix 
elements. In general, as noted before, relaxation also re- 
sults in dephasing, and in particular for three-level ladder 
systems such as the one discussed in this experiment both 
Fio and r2i contribute to the the off-diagonal decay rate 
721 = 712 ■ 

Using now the master equation p = —{i/H)[H^'^^\p] + 
C'^^^^Ip], with the effective Hamiltonian Eq. ((56)) and the 
effective Liouvillean Eq. ([55)) we are able to provide a 
complete characterization of the dynamics of the system. 



IV. SPURIOUS EXCITATIONS TO HIGHER 
LEVELS: A 5-LEVEL MODEL 

In the previous sections we have developed a three-level 
model with dissipation for our phase qubit. Truncating 
the Hilbert space to three levels is a good approxima- 
tion, as one can see when comparing the predictions of 
this model with the experimental data (next sections) . In 
addition, as we have shown in the previous section, ana- 
lytical results can be derived. However, due to the fact 
that the phase qubit is a multi-level system with finite an- 
harmonicity, one source of errors for the three-level model 
comes from leakage to the higher excited levels. This is 



phasing part is unchanged in the doubly- rotating frame. 
Thus decoherence is described by the following effective 
Liouvillean, 

J6{1,2} 

7^^ 

- E -^'^jjP'^kk, (28) 

j,ke{0,l,2}:j^k 

or in explicit matrix form 



(29) 

I 

due to the off-resonant coupling of the fields with the 
higher-level transitions and subsequent transitions into 
the first three levels. This effect can be accounted for by 
a straightforward generalization of the three-level model 
to more levels, and then numerically solving the corre- 
sponding master equation. For our simulations we chose 
five levels. We found that further increasing the number 
of levels produces no significant change in the results. To 
avoid cumbersome equations, we only briefly describe the 
model here. We follow exactly the same procedure de- 
scribed in the previous two sections: starting from Eqs. 
(|T2)) and Eqs. ()13p . we calculate the matrix elements of 
the perturbed harmonic oscillator H' = Ho + in 
the basis of Hq- The resulting 5x5 matrix is diagonal- 
ized numerically; the time-dependent term A$$i-f is then 
added and the full five-level time-dependent Hamiltonian 
H{t) corresponding to Eq. (fT7)) is obtained. For the dis- 
sipation, we use the straightforward generalization of the 
previous Liouvillean, 

^Ps] ^ E -^Y^i'^'^J-hjPscrjJ^l-crjjPs- PsCTjj) 
je{l,2,3,4} 

7'^ 

J2 ^o-jjPsc^fcfc, (31) 

j,ke{0,l,2,3A};j¥'k 

Note that in Eq. (|3T)) we have neglected terms that mix 
the interlevel relaxation rates - of the type appearing as 
the last one in Eq. (|55)) . The justification for this is 
that, similar to the approximation involved when going 
from Eq. ^7}i to Eq. ([28]) in Section [HTI these terms 
will average out when moving to a multi-rotating frame 
(defined by extending Eq. (|M| to five levels). Thus we 
can neglect them already in the Schrodinger picture. 
Finally, we put the master equation {d/dt)ps(t) = 



ZFioPii 

-(rio + 7fo)p 
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-(Fio 4 
-SFiopii 



lw)Poi 

+ 2r2lP22 



-(r2i +72o)P20 -(Fio + r2i +7|i)p: 



~(r2i - 
^(rio + r2i 



l2o)PO 



l2l)P 
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— (i/H) [H{t), ps{t)] +C[psit)] in matrix form and solve it 
numerically. 

V. DETERMINATION OF THE DEPHASING 
RATES 

We now address the question of how to determine the 
numerical values of the parameters that enter in the 
Hamiltonian and in the Liouvillean of the models pre- 
sented. First, a number of independent experiments are 
performed that allow us to find ajio,W2i,rio, and r2i. 
The transition frequencies wio = Stt x 8.125 GHz and 
W21 = 27r X 7.975 GHz, are determined immediately from 
standard single-tone and two-tone spectroscopy (see Fig. 
13]). The relaxation rates are determined as follows: for 
Fio we apply a tt pulse resonant to the first transition, 
taking the system from |0) to Then we let the qubit 
decay and determine the occupation probability by delay- 
ing the measurement pulse with respect to the tt pulse. 




where we use the same notations as before, Ap = wiq— Wp, 
and 710 = Fio -t- ^fg. In the steady state p = 0, and 
the occupation probability of level as a function of 
detuning, is obtained as 

4A^rio + 710 i^^p + rio7ioJ 

This probability reaches its maximum value at ujp = loiq. 
The width SfiQ — Sujio/27r of this spectroscopy peak. 

Slow = 27rJ/io = y^7?o + 2f^^7io/rio « 710, (34) 

is determined by solving the equation Pi(Ap = ttJ/iq) — 
Pi(Ap = 0)/2. The result shows, as expected, that for 
relatively small probe fields the width of a spectroscopy 
line of a two-level system is given by the sum of the re- 
laxation rate and the pure dephasing rate. Using Eq. 
([M)) we can now determine jf^. We use the values of 
the independently-measured decay rate Fio = 27r x 7 
MHz and find SfiQ = 14 MHz from the black contin- 
uous line (the measured data) in Fig. [31 This yields 
7fo w Slow - Fiq = 27r x 7 MHz. 

To obtain the other dephasing rates J20 ^21 there 
is no simple analytical result that we can use, there- 
fore we need to rely on simulations and direct fittings 
with spectroscopy traces involving level |2). We write 
J21 = 7fo + 720 ~ ^' ^^'^ ^^'^ thai the best fittings are 
obtained for ^20 = 27r x 16 MHz and e = 27r x 5 MHz 
{•^21 = 27r X 18 MHz). As we show in Appendix iBl e is a 



This probability decays exponentially as a function of 
time, allowing us to extract the relaxation Fio = 27r x 7 
MHz. The same technique can be applied to record the 
decay of the |2) state to the |1) state, while simulta- 
neously driving the first transition continuously to en- 
sure a non-zero population on the state |1), and we get 
F21 = 27r X 11 MHz. 

The pure dephasing rates 7fg, 7^0, and 7^^ can be ob- 
tained from the spectroscopic linewidths in the following 
way. We first do a single-tone spectroscopy, measuring 
the sum of the occupation probabilities of the excited 
states at zero coupling field (black solid curve in Fig. [3]) . 
Since Jlc = in the first approximation we may treat the 
system as a quasi-two-level system under monochromatic 
driving, and neglect the leakage to higher levels. This 
experiment thus represents a simple spectroscopy of the 
two-level system {|0), |1)}. We, therefore, can write the 
master equation with the rotating-frame effective Hamil- 
tonian Eq. (|26)) and Liouvillean Eq. ((29)) truncated to 
the ground state and first excited state and obtain 



(32) 

I 

0.1 6 1 1 1 1 1 1 1 




(Op(coJ/27t(GHz)' 

FIG. 3: (color online). The spectroscopy traces (solid curves) 
taken from Fig. 1(c) in Ref. [ij. The (black) solid trace is 
the single-tone spectroscopy experiment on the effective two- 
level system showing the |0) — >■ |1) spectral line centered at 
8.135 GHz. The (red) solid trace is a two-tone spectroscopy, 
revealing the |1) — >■ |2) spectral line centered at 7.975 GHz. 
The dotted curves are fits to simulations of the full five-level 
master equation. The value of the probe field Rabi frequency 
used was fip = 27r x 3.5 MHz. 

measure of correlation between the fluctuations of level 
1 and level 2. Reasonably good fittings are obtained also 
for e = MHz (see Fig. 1^). As we will see in the 
next section, the advantage of the approximation £ « 
MHz is that relatively simple analytical results can be 
obtained. 

In Table U we summarize our results for the parame- 



ii^pipoi - Pio) + 2TiQpii ~i[^lp{pii - Poo) - 2ApPoi] - 7ioPoi 

-i[Op(poo - pii) + 2Ap/9io] - 710P10 -i^pipoi - Pio) - 2riopii 
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ters of the three-level model. Even better fittings can be 
obtained if one employs the five-level model described in 
Section llVi thus accounting for the leakage outside the 
subspace spanned by {|0), |1), |2)}. For this subspace we 
now have all the numerical values needed to simulate the 
dynamics. Now, in the Hamiltonian, the additional ma- 
trix elements corresponding to the fourth and fifth states 
(1 3) and |4)) can be calculated based on the information 
already obtained from the three-level model. The only 
unknowns are the dissipation parameters associated with 
the last two levels, for which an independent experimen- 
tal determination is not easily available. A rough esti- 
mation of these parameters is based on the result that, 
since dissipation is due to coupling to the external elec- 
tromagnetic environment, we can expect the same scal- 
ing to continue. This estimation does not need to be too 
precise: we have checked numerically that even relatively 
large errors in the estimation of these higher-level dissi- 
pation parameters do not result in significant mismatches 
with the experimental data. The results of simulating the 
five-level model are plotted with dotted lines in Fig. [31 
The dotted red curve fits almost perfectly the |0) O |1) 
transition, and a very good fit (black dotted lines) is ob- 
tained for the |1) O |2) transition. Thus, we conclude 
that the five-level model confirms the predictions of the 
three-level model and further improves the fitting with 
the experimental data. 

VI. ANALYTICAL RESULTS FOR THREE 
LEVEL SYSTEMS 

In this section we show that relatively simple analyti- 
cal results can be obtained if we use the approximation 



721 « 710 +720, where 710 = rio+7u)' and 720 = r2i-l-7|'o 
are the total off-diagonal decay rates introduced before 
in Eq. ((30)) . This relation is satisfied reasonably well for 
our system, as we show in Appendix [Bj Note that al- 
though there is no a priori reason for the general validity 
of this equation in a three-level system, we do expect it 
to be a good approximation for many cases of interest. 
Indeed, the error that we make by this approximation in 
the total off-diagonal decay rate 721 is e = 7io+72o~72i, 
which should be compared with 710 -I- 720- But the lat- 
ter quantity is relatively large, since it contains the sum 
of two relaxation rates and two pure dephasing rates, 
710 + 720 = r2i -I- Fio -I- 7^) + 7fo- For example, in our 
case we have 710 + 720 = 27r x 41 MHz and e = 27r x 5 
MHz, meaning that our approximation results in an over- 
estimation of about 13% of the true value of 721. 

Using Eqs. ([^5)) and ([^5)) we can now solve for the 
density matrix p, which satisfies p = —{i/h) [iJ^^^^p] -|- 
£('=ff)[/9]. In the steady state, = 0, and the ele- 
ments of the density matrix can be found in analytical 
form using Mathematica. The full solution can be put in 
a simpler form if some further approximations (see also 
Appendix [C]). corresponding to the actual values used in 
the experiment, are used. 

The first approximation that we employ is that the 
ratio between fl^ and any product of two 7's is much 
smaller than 1 {Q^ <C 720, 7iO' 7io72o); the largest ratio 
for our experiment is ^I^/^iq = 0.06. Next, we neglect 
the square of the probe field with respect to the square of 
the coupling field {flp/il'^ <^ 1); for the parameters used 
in our experiment, this ratio is smaller than 0.02. 

Under these conditions, we obtain the stationary oc- 
cupation probability of level |1), 



nl [4A2710 + 720 {nl + 710720)] 
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16A4 + + 710720)' + 4A2 (72^ + 7I, - 2n^J 



(35) 



This analytical equation agrees well with the exper- 
imental data and also with a more complete five-level 
model described previously. For Ap = 0, i.e. at the 
frequency where in the absence of the coupling field 
there was an absorption peak, Eq. (|35|) yields Pi = 
l20^io^pi^l+liol2o)^^ , showing that the population of 
level 1 decreases to zero as ftc increases. As a function of 



Ap (or Ljp) the occupation probability has a double-peak 
shape (see later in Fig. ([S]). We will refer to these peaks 
as the Autler-Townes peaks. A more detailed analysis is 
presented in Section IVlII 

The splitting Sf at = 6ujat/'2t^ is twice the value of the 
positive solution of the equation d^pPi = 0; we get 



. xf /fic 710 + 720 +710720 - 720 +710720) 

dujAT = 27r(5/AT = \ • (36) 

V 710 



When comparing this formula with the experimental data the agreement is very good (see Section lVHi Fig. O . One 
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parameter 


[Ji,i_i (GHz) 


r,,,_i (MHz) 


Jfo (MHz) 


i = 1 


27r X 8.135 


27r X 7 


27r X 7 


i = 2 


27r X 7.975 


2-K X 11 


2tv X 16 



TABLE I: Parameters for the three-level model. For 'yl^-^ we 
find = 27r X 18 MHz. 



also notes that if ilc is much larger than the decoherence 
rates the dependence is linear Sloat ~ ^c', otherwise, dis- 
sipation tends to move the Antler- Townes splitting away 
from linearity. 

The matrix elements poi = Pio can be calculated under 
the same approximation of a weak probe field as above, 
with the result 



Pw 



fip(2Ap - z72o) 



-4A2 + 710720 + 2iAj,(7io + 720) ' 



(37) 



Although this off-diagonal matrix element is not directly 
measurable in our experiment, it plays an important role 
in other experimental configurationsi^i^^. 

The rest of the elements of the density matrix are 
too complicated to be listed here. Still, we can make 
progress by restricting ourselves to the case of coupling 
fields larger then all the linewidths (more precisely, the 
case in which the products 7|g, 7^0 and 710720 can be 
neglected with respect to fie)- The expressions for the 
elements of the density matrix in this case are given in 
Appendix O 



VII. COMPARISON WITH THE 
AUTLER-TOWNES EXPERIMENT 

In the previous sections we have described how to ob- 
tain the parameters entering in the Hamiltonian and in 
the Liouvillean, which determine the dynamics of the 
qubit (for the three- level model, see Table HI. In this 
section we describe how to use this information for mod- 
eling the behavior of the system at large coupling pow- 
ers, which produces the Autler- Townes effect. This effect 
is seen in the spectroscopy trace as the splitting (some- 
times called dynamic Stark splitting) of the line corre- 
sponding to the first transition into two peaks. In the 
doubly-rotating frame, this can be understood as a con- 
sequence of the dressing of the transition |1) — >■ |2) by 
the electromagnetic field ilc to which it is coupled (Fig. 
H]). Effectively, the weak probe field fip sees the doublet 
of the two dressed states formed around the qubit state 
The dressed-state picture follows naturally when the 
fields are treated quantum-mechanically^, but a similar 
picture emerges also in our formalism, which uses classi- 
cal fields. From Eq. ((28)) . by neglecting the probe field 
and with the coupling field at resonance uJc — t^2i, the 
two nonzero eigenvalues form a doublet separated by an 
energy hflc (which is Sujat at Qp — 0). 

In Fig. [5]we show the same data as in Fig. 3(a) in Ref. 
prj . The (black) solid curve corresponds to fie = 27r x 36 
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FIG. 4: Schematic of the Autler- Townes effect as a conse- 
quence of the formation of dressed states for the |2) o |1) 
transition due to a resonant couphng field. 



MHz and the red one to ftc ~ 2Tr x 66 MHz; the value of 
ilp was 27r x 3 MHz. We first fit the experimental data 
by using the analytical model with three levels presented 
in Section mil fsee Eg. (j^S))). The results are shown with 
(green) dashed lines in Fig. [S) Then we use the nu- 
merical solution of the full five-level model to obtain the 
dotted lines in Fig. [SJ and one sees that the fit is almost 
perfect. When plotting all the numerical results one has 
to remember to add the residual tunneling amplitude of 
approximately 0.057, which was mentioned in Ref. jT]| . 
and which is also visible in the off-resonant part of the 
black trace in Fig. |21 This is specific to the measuring 
scheme (based on switching of a nearby SQUID). This 
offset also has to be added when comparing Eq. ((35)) 
with the measured data. At high coupling powers, the 
field flc will couple to the |0) — > |1) transition as well, an 
effect that is not captured by Eq. ([55)) . As a result, when 
using Eq. ([35)). we expect that at low powers the value 
of this offset is close to 0.057, but at high powers (e.g. 
ilc = 2tt X 66 MHz) the offset value has to be adjusted 
appropriately (see Fig. 

Another observation regarding the three-level model 
is that it predicts symmetric Autler- Townes peaks, while 
the experiment and the full model give peaks with slightly 
different heights and displaced to higher frequencies (see 
Fig. [S|). However, the distance between these peaks (the 
Autler- Townes splitting) is predicted remarkably well in 
this simple model. Indeed, in Fig. [6|we show a compar- 
ison between experimental data (as in Fig. 3(c) in Ref. 
ll| ) and the results for splitting as given by Eq. (|5S)) . 
We notice that the three-level model with dissipation pre- 
dicts a bending of the Autler- Townes frequency splitting 
at low powers, which is observed in the experiment and 
is not captured by the analysis of the energy levels in the 
three-level model without dissipation. 



Finally, we can simulate the case when both the cou- 
pling and probe frequencies are swept, and we get again 
a very good agreement with the data, as in Fig. 4 of Ref. 
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FIG. 5: (color online). Autler-Townes splitting in a phase 
qubit. The experimental spectroscopy traces (solid lines) are 
taken from Fig. 3(a) in Ref. 11]. The dotted lines are simu- 
lations of the master equation corresponding to the full five- 
level model, with the parameters previously extracted from 
the spectroscopy data. The (green) dashed lines are plotted 
by using Eq. p5|) with the same values (the offsets of base 
level for flc ^ 2tv x 36 MHz and Qc = 2tt x 66 MHz are 0.057 
and 0.105, respectively). Note that the inclusion of higher 
levels at large values of the coupling results in asymmetric 
Autler-Townes peaks and in the displacement of the spec- 
troscopy curve toward higher frequencies. The vertical (blue) 
dashed line indicates the first transition uio = 2n x 8.135 
GHz. At this frequency the lower curve (taken at Qc ~ 27rx36 
MHz) has a minimum, while for the higher curve (taken at 
Qc = 2-kx 66 MHz) the minimum, marked by the blue vertical 
arrow, is displaced by 10 MHz to a value of 8.145 GHz. 



VIII. CREATION OF DARK STATES 

An important fundamental issue is whether true dark 
states and electromagnetically induced transparency can 
be created in systems of superconducting qubits. A dark 
state is the zero-eigenvalue eigenvector of the effective 
Hamiltonian H^'^^^ of Eq. (|28|l for the resonant case, 
Ap = Ac = 0. The Hamiltonian is diagonalized as 

|D) = cose|0) - sine|2) eigenval. =0, 
^ [sine|0) ± |1) + cose|2)] eigenval. = ± ^O^ + ni, 

where the mixing angle Q is given by tan0 — ilp/ilc- 
The first state is characterized by the absence of pop- 
ulation on |1); contrary to this, the other two eigenval- 
ues have a considerably larger population of the state 
|1). Moreover, in the limit Vlp <^ flc (which corresponds 
well to our experimental parameters), we notice that the 
dark state becomes approximately equal to the ground 
state |0), and the other two eigenvectors become approx- 
imately (l/-\/2)(±|l) -I- |2)). Thus, ground-state ampli- 
tudes are 1 for the dark state and for the other two 
eigenstates. These features are a hint that the dark state 
is easily distinguishable from the other two eigenstates. 




n^/27t (MHz) 



FIG. 6: (color online). The Autler-Townes frequency split- 
ting SujAT I (Stt) as given by the simplified three-level model 
with dissipation Eq. p6|) . shown as a (red) dotted line. The 
straight line is the prediction of the three-level model with- 
out dissipation. The experimental data (circles) are from Fig. 
3(c) of Ref. iTl[. 

Indeed, by looking at the experimental data presented 
above (small populations for the levels 1 and 2) it is clear 
that they are consistent with the stationary state of the 
system being close to a dark state. In order to quantify 
this observation, we calculate the distance between the 
steady-state p*^^*) and the dark state \D) by the fidelitjs^ 

J-p)[p(^*)] = 7(D|p(^')|D) (38) 

cos 26 (st) (st)N sin 29 (st) , (st)^ 
= (Poo -P22 ) —\PiQ +P\n ) 

+ (39) 

using fip and fie as adjustable parameters. The results 
are presented in Fig. d?]). We have also checked that the 
same result is obtained in the Schrodinger picture (the 
non-rotating frame) by evolving numerically the system 
for a time much longer than all the timescales set by dis- 
sipation and driving, and calculating the fidelity with re- 
spect to the dark state in the non-rotating frame, where 
it takes the form cos0|O) + exp[— j(ajc + Wp)t] sin 0|2). 

Similarly, one can calculate the state purity Tr{(p(''*)) }; 
the result is plotted in Fig. [51 As expected, for rela- 
tively low amplitudes of the probe pulse, as used in our 
experiments, the state is very close to a pure state. 

For our experiment the fidelity reached values above 
0.95 (see Fig. [T]) at pumping fields fie higher than 27r x 50 
MHz. Thus, to a reasonable approximation, these exper- 
iments can create steady states that are close to ideal 
dark states. Steady states could be useful in future 
quantum processors for state preparation and as sources 
of entanglement^. A more in-depth characterization of 
these states, and the effect of relaxation and dcphasing, 
is given in Appendix[Bj One can anticipate, for example, 
that the decay from the level |2) will inevitably decrease 



11 



Fidelity of the steady state 
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FIG. 7: The fidelity of the steady-state as a function of the 
pump and probe amplitudes. The probe and field amplitudes 
are directly related to the respective Rabi frequencies by Q,p — 
0.69gp and fic = <?c. The blue line corresponds to Qp — 2tt x 
3.5 MHz {qp = 2n x5 MHz), and to the range of experimental 
values for 0,^ presented in Fig. [6l 



Purity of the steady state 
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FIG. 8: The purity of the steady-state as a function of the 
pump and probe amplitudes. As in Fig. [71 the blue line 
corresponds to fip = 27r x 3.5 MHz {gp — 2n x 5 MHz) and 
a range of experimental values for Qc as presented in Fig. [S] 
Also, the probe and field amplitudes are directly related to 
the respective Rabi frequencies by Q,p = 0.69(?p and fie = gc- 

its population, and will "contaminate" the true dark 
state with an incoherent population on level |1). This ob- 
servation is confirmed by the analysis in Appendix [CJ A 
somewhat more advantageous situation could be achieved 
if the system is used in a A-configuration, in which case 
the lifetime of the upper-state component of the dark 
state is increasedi^. Yet, even for the ladder configura- 
tion discussed here, the stationary state is close to an 
ideal dark state. 

In general, the main obstacles in increasing the fidelity 
are the relatively short (compared to atomic physics) de- 



coherence times. However, the advantage offered by ar- 
tificial systems is that larger couplings with the fields 
can be realized: this allows experimentalists to increase 
the strength of the field until the Autler-Townes splitting 
well exceeds the natural spectral linewidth. 



IX. CONCLUSIONS 

We have studied a phase qubit under a two-tone mi- 
crowave irradiation applied close to the |0) — )■ |I) and 
|1) 1 2) transition frequencies. We provided a de- 
tailed quantum-mechanical theoretical description of the 
three-level system interacting with the fields, including 
relaxation and dephasing. The relaxation and dephas- 
ing rates were extracted from independent experiments. 
These parameters were then used in the analysis of the 
Autler-Townes effect to compare the measurements with 
the theoretical predictions. We found that the value of 
the Autler-Townes splitting and the occupation probabil- 
ities predicted by the three-level model are in very good 
agreement with the experimental data. This agreement 
improves if the next higher levels are included and full 
numerical simulations are performed. Finally, we ana- 
lyzed the structure of the stationary states thus created, 
and show that they are very close to true dark states. 



Acknowledgments 

G.S.P. thanks Petr Anisimov for useful discussions. We 
acknowledge financial support from the Academy of Fin- 
land (no. 129896, II8I22, 130058, and 135135), from 
the National Graduate School of Material Physics, from 
NIST, and from the European Research Council (StG). 



Appendix A: Energy relaxation 

We define the density matrix of our whole system 
{system S ® reservoir R) as xi^)- Since our interest 
centers only around the resulting state of the system S, 
irrespective of the outcome of R, we want to obtain an 
equation of the reduced density matrix p{t) = Trulxit)] 
in order to describe the dynamics of the system S induced 
by the evolution of the total system. The aim is to derive 
a so-called Markovian master equation for the reduced 
density matrix p following CarmichaelSS,, and Breuer and 
Petruccione^^. 

The Hamiltonian of the whole system is assumed to 
have the general form 

Htot ^Hs + Hr + Hsr, (A1) 

where 

Hs = hujw\l){l\ + h{ujw + ^2i) |2)(2| (A2) 
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denotes the Hamiltonian of a three-level system in its en- 
ergy eigenbasis {|0), |1), |2)}, with transition frequencies 
wio, and uj2i. The Hamiltonian Hn^ 

HR = h^uJaaiaa (A3) 

a 

is a collection of harmonic oscillators with frequency uja 
and corresponding creation (annihilation) operator aj^, 
(oq,), and 



HsR = K (/i"|0)(l| + /2a|l)(2|) + h.. 



(A4) 



describes the coupling between the inter-level transitions 
of the three-level system and the reservoir oscillators. 
The coupling constants fja {j = 1, 2) are assumed to be 
small (weak-coupling limit). Note that we do not take 
the |0) ^ 1 2) transition into account due to the small 
anharmonicity of the system. 

We transform the Schrodinger equation for the density 
operator of the whole system x, (from now on we take 
H=l) 



X = ~i[Htot,x], 



(A5) 



into the interaction picture, with x(t) — exp[i{Hs 
Hji)t]x{t) exp[—i{Hs + Hii)t], and obtain 



X = -i[HsR{t),x], 



(A6) 



where HsR{t) is given by HsR{t) 
HR)t]HsRexp[^t{Hs + HR)t]. 

By integrating Eq. (|A6p and substituting for x(t) in- 
side the commutator in equation (jA6p , we get a Liouville 
equation in the integro-differential form 



-i[HsRit),xm 



which assumes that the influence of the system on the 
reservoir is negligibly small due to the weak coupling be- 
tween them. Therefore, the total density matrix x{t) at 
time t can be expressed as a tensor product 

X(i) -/5(<)®i?(0), (AlO) 
with i?(0) = cxy>[-I3Hr)/Z, where fi ^ l/fcsT, and the 
partition function is Z — Tr[exp(— /3i7jj)]. 

Hence, tracing both sides of Eq. (jA7p over the reservoir 
variables, and eliminating the term —iTTR[HsR{t),x{^)] 
with the assumption Tr r[H s R{t) R{0)] = 0, we write Eq. 
(|T7)) as 

m = -f dt'TvR { [ilsRit), [HsRit'), Pit') ® i?(0)]] } . 

(All) 

In general, the evolution of the system depends on its 
past history (p(t') in the integral). However, it is known 
that for a reservoir with many degrees of freedom this 
memory can be erased [Markov approximation), which 
allows us to make the replacement 

Pit') ^ Pit), (A12) 
and leads us the following master equation 

m = - y dt'TvR { [ilsRit), [HsRit'),pit) (g> RiO)]\ } . 

(A13) 

To get a more explicit form for the master equation, 
we need to insert HsR^ into Eq. (jAlSp . so we rewrite Eq. 
^ as 



H 



SR 



h.c. 



(A14) 



- / dt' 
Jo 



HsRit),[HsRit'),xit')] . (A7) where aj.k = (Jjk = are operators for the sys- 

tem and r| = fjaO-a operators for the reser- 



The reduced density matrix for the system S in the in- 
teraction picture is given by 



Pit) = TrRxit) = e'"''^*pit)e 



-iHst 



(A8) 



voir. Transforming into the interaction picture by 
(Tj_ij(t) = expiiHst)aj-ij expi—iHst) and f](i) — 
exp(ii/fli)r^ expi—iHRt), all these operators become 



We assume that the interaction between S and R is 
switched on at time t = 0. So at this initial time, S 
and R are uncorrelated and the total density matrix is 
factorized as the direct product 



and 



X(0) = p(0)®i?(0)=x(0), 



(A9) 



where i?(0) denotes the initial reduced density matrix for 
the reservoir. 

In order to eliminate x from Eq. (jA7p . we perform 
the first approximation, known as Born approximation, 



The resulting master equation in the Born-Markov ap- 
proximation is then 
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P = 



E 
E 



df 



{t)m) + h. 



c, o-fc-i 



fc(Ofl(t') + /i.c.,p( 



i?(0)]]} 



+ [dk,k-i{t')pdj,j-i{t) - dj,j^i{t)dkM-i{t')p] {fj{t)fk{t'))R 
+ [dj^j-i{t)pdk,k-i{t') - pdk,k-i{t')aj,,-i{t)] {Tk{t')tj{t))R 
+ [(Tk,k-i{t')pd--i-ij{t) - aj-ij{t)ak,k-i{t')p] (f](t)ffe(t'))i? 
+ [aj_ij(t)pCTfc,fc-i(t') - pa-fc,fc_i(t')o-j-ij(i)] (ffe(t')f](t))i? 



(i)pafc_i,fc(i') - pak-iAt>j,i-i{t)] (fl.(i')f, (^)>i^} , 



r 



(A15) 



where 





= Trfl 


R{Q)r]{t)rl{t') 




= Trfl 


'Rio)rlit')f]it) 


f,Wf,(t'))i^ 


= Trn 


'R{0)r,{t)Tk{t') 


ffc(t')f,W)« 


= Tr;^ 


"i?(0)ffe(t')f,W 


f](t)f,(t'))« 


= Tr^ 


'Rio)r]{t)fkit') 


ffe(<')f]W)« 


= TrR 


'Rio)rk{t')f]it) 


f,(t)fi(t'))« 


= TvR 


'Rio)r,{t)flit') 


fI(t')f,W)fl 


= TrR 


"^(0)fl(t')f,(i) 



are two-time correlation functions. These two-time cor- 
relation functions are expected to vanish for t — t' larger 
than a time scale Tc oc h/ksT called the correlation time 
of the reservoir. Thus, the integral in Eq. (jA15l) is 
nonzero only for times < t — t' < Tc, and for t' outside 
this time interval, p(t') should not affect p{t) at time t. 
As discussed in Ref. \2&\ we assume that the time scale 
for significant change of the system is much larger than 



the reservoir correlation time; then the Markov approxi- 
mation Eq. (|A12I) holds. 

Since we also suppose that the reservoir is in ther- 
mal equilibrium, the correlation functions with terms 
{'^a'^a)R '^^ {<^aaa)R vanish, and the terms {a]^aa)R and 
(aQflDfl give 

TvR[R{Q)aia^] = n{Lo^,T), (A16) 
TvR[R{Q)a^ai] = l + nK,T). (A17) 

Here n((u;c,r) has the meaning of mean photon number 
for an oscillator with frequency Ua in thermal equilibrium 
at temperature T, 



1 



(A18) 



We change the summation over the reservoir levels in 
the nonvanishing correlation functions to integrations by 
introducing a density of states G{uj), and use t = t — t' 
as a new variable in Eq. (jA15l) . By substituting the 
nonvanishing correlation functions and their correspond- 
ing system operators in equation (lAlSp . the Markovian 
master equation in the interaction picture is explicitly 
written as follows (define the detuning A = wio — W21) 



'p = (aiopaoi - (Toiaiop) / dr / dwe'^^-'^'°^^G{uj)\h{uj)\''n{<^,T) 

Jo Jo 



+ (^oiMo-Moaoi) / dT dwe^('^-"-)^GH|/iHnn(c^,r) + l] 
^0 Jo 

pt poo 

+ (^2ipaoi - aoiaaip) / dr duje'^'^-'^''^^e-''^'G{u;)fi{uj)f;{u;)n{uj,T) 



^0 
t 



+ (aoipa2i - /5a2iaoi) / dr dwe*('^-'^-)^e-^^*GM/iM/2*M[n(c., T) + 1] 
Jo Jo 
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+ KiMo - (Tioaoip) / dr dc^e-*("-"-)^GH|/iHnn(c^,r) + l 

^0 ^0 

ft pOO 



+ (^TioP'Toi - paoiaio) / dT / duje-'^^-^'°'^^G{Lo)\h{u:)\''n{Lo,T) 
Jo Jo 

/•i poo 



Jo JO 

+ (^710^2 - ^Ti2aiop) / dr /" dwe*('^-'^-)^e^^*GH/r H/2Hn(c^, T) 

pt pOO 

+ {<Ji2P<Jio - P<yio<Ji2) dT dwe*('^-'^-)^e^^*G(L.)/rM/2H[n(w,r) + i] 
Jo Jo 

pt poo 

+ {<T2ipai2-<Ti2<T2i~p) / dr / du:e^^-^'^>G{uj)\f2{u:)\''n{Lo,T) 
Jo Jo 

pt poo 

+ {<Ji2P<J2i - P<J2i<Ji2) dr dwe*('^-'^")"GH|/2HHn(c^,r) + l] 
Jo Jo 

pt pOO 

+ {<Jmpi2i - cT2ia„ip) / dT / duje-'^'^-'^^°^^e-'^'G{u)h{uj)f;{u)[n{u,T) + l] 
Jo Jo 

pt pOO 

+ (^2ipaoi-/5aoia2i) / dr dioe~''-^~^-°')^e''^'G{oj)h{u:)f;{uj)n{u:,T) 
Jo Jo 

pt poo 

+ {<Ji2P'J2i - <T2i'Ji2p) dT dc^e-^('^-'^-)^GH|/2Mnn(c^,T) + l] 



Jo 

t poo 



+ {<T2iP<Ji2-p<Ji2<J2i) dT dwe-*('^-'"^i)^GH|/2Hpn(cj,T). (A19) 
Jo Jo 



I 

Since t is the typical time scale for changes in p, which gral to infinity, and then we evaluate the integrals in Eq. 
is much longer than reservoir correlation time Tc, we can (IA19P as follows (without loss of generality, assume /i 
approximately extend the upper limit of the time inte- and /2 are real) 



t /"OO 

dT / dwe±*('^-'^-)^G(L.)/f (w)n(L., T) = ^G(c^io)/? (c^io)^(^io, T), 
Jo 

t poo 

dT / dwe±*('^-"-)^G(c^)/2'(w)n(c^,T) = ^G(c.2i)/2 ('^2i)fi(^2i, T), 
Jo 

t /"OO 

dT / dwe±^('^-"-)^G(^)/i(L.)/2(w)n(c^,T) = nG{ujio)fi{^io)f2{uJio)n{oj,o,T), 
Jo 

t POO 

dT / dwe±*('^-'^-)^G(c^)/i(c^)/2(c<.)n(c^,T) = ^G(c^2i)/i(c^2i)/2(c^2i)n(L.2i,T). (A20) 
Jo 



I 

Since ojiq ~ UJ21 (wio = 27r x 8.135 GHz and aj2i = the master equation in following form 
27r X 7.975 GHz), we can make the approximations 

G(wio) ~ G(w2l) = G, /l(2)(wio) ~ fn2){^2l) = /l(2), 

and n(tJio,T) « n(w2i,7') = ■ Tio., 

p = —n{2aiQp(Joi - aoop - ptroo) 

Combining equation (|A19p with equation (|A20|) . and Pio 
using aoia2i = cr2icroi = criocri2 = cri2crio = 0, we get + —{n + l)i2aoipaio - (Jiip - pail] 
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—n[2a2ip(Ti2 ~ (Jiip - pail) 



(A21) 



21 



(n + l)(2cri2po-21 - 0-22/5 - /5cr22) 



Ke ''^*croi/5cr2i + *'^*(ctoiPo-21 + cr2ipo-oi) 

j Af ~ 

Ke cri2/9Crio + K1 



^*(o'i2/5fTio + crio/5cri2), 
2^G/2, T21 



where the energy relaxation rates Fio 
2TTGfl and k = V^w^^- 

Transforming back to the Schrodinger picture, and tak- 
ing the zero temperature hmit n — >■ 0, we finally obtain 
the Markovian master equation 



in the final form [see also Eq. ((22|) ] 

P = -iiHs, p] + k(o-01/5(T21 + (Tl2/5Crio) 

— ^(2o'oipoio - diip - peril) 

H ^(20'12P0-21 - 0-22P - PO-22). 



(A22) 



(A23) 



In a doubly-rotating frame, terms with k average out to 
zero on the time scale of our measurements, due to the 
factor e^^''* (see Sec. Hill) . 



Appendix B: Pure dephasing 

In this Appendix we derive the dephasing part of the 
Liouvillean for a three-level system, as introduced in Sec- 
tion [TTTl We also present numerical results bearing on the 
validity of the approximation 721 ~ 720 + 7io- For our 
system we suggest that an appropriate model for dephas- 
ing could be that of a harmonic oscillator with fluctuat- 
ing frequencies in addition to uncorrelated noise due to 
(energy-conserving) virtual transitions with the environ- 
ment. 

For two-level quantum systems, it is known that de- 
phasing is caused by longitudinal noises (i.e. along CTz), 



J 



which induce fluctuations in the transition frequency. Al- 
though for three-level systems we do not have a Bloch 
sphere picture to help our intuition, we can still develop 
the theory by using the matrices ajk defined in Appendix 

m 

We consider the ladder three-level Hamiltonian Eq. 
^H), and 



Hfi{t) = nSuio{t)\l){l\+n[5uj2i{t) + Suio{m2){2\ (Bl) 

describing random fluctuations of transition frequencies. 
The total Hamiltonian in matrix form reads 



ujio + Sujioit) (B2) 

UJ2I + i^lO + ^^2o{t) 

where we define 6uj2Q{t) — Si02i{t) + SuJio{t). 

In the interaction picture, the dynamics of the three- 
level system's density operator, 



p{t) = cxp{iHst)p{t)exp{-iHst), 



(B3) 



follows from the Schrodinger equation and a perturbation 
expansion to the second order in noises. 



ih'p = [Hfi{t),p{t) 

« \Hfi{t),p{0) 

Hfi{t),m 
dt' 



1 

ih 



dt' 



Hfi{t'),p{t') 



Hfiit), Hfi{t'),p{t' 



(B4) 



where Hfi{t) — exp{iHst)Hfi{t) exp{—iHst). The dou- 



ble commutator 



from the last 



[Hfiit), [H flit'), Pit') _ 
line of Eq. (|B4I) can be calculated and has the follow- 
ing matrix form 



poi{t')[SLJlo{t)SLUlo{t')] Pa2{t')[SLU20{t)SuJ20{t')] 

pioit')[6wio{t)6ujio{t')] pi2{t')[Sujio{t) - SLJ2omSujio{t') - dLJ2o{t')] 

P20it')[5uj20{t)6LO20it')] P2l{t')[SuJioit) - 5^20 (0] ['5^10 (*') ^ SLU2o{t')] 

m 



Then, after averaging over the fluctuations, Eq. (jB4p can be rewritten in the following matrix form, 
I 



Hfiit),m 



16 





-7foPio(i) 





Sl2 - S2l)p2l{t) 



10 



720 



71^0^02(0 

S'l2 - S2l)pi2{t) 





(B6) 



where we have as usual the density matrix element as 
Pjk(t) = 0|y5(t)|fc) and we made the approximation 
p{t') « p{t) (see also Appendix IX)) ). The dephasing rates 
are defined as 



l2Q 



dt'{Sujio{t)6ujio{t')), 



dt'{Su;2oit)Suj2oit')), 



and the cross spectral densities as 



Si2^ J dt' {Sujio{t)5uj2o{t')} = S'21. 



(B7) 
(B8) 

(B9) 



Let us now define define a parameter e = S12 + S21, which 
parametrizes the effect of correlations between levels 1 
and 2. We can now write the dephasing terms in Eq. 
(lB6l) as 



+ 2 (^llP'''22 + a'22pcrii) , 



(BIO) 



or equivalently in the form used in the main part of the 
paper Eq. (^5]) . 



^dep[p] 



E7ifc 



(JjjPCTkk, 



(Bll) 



j,fe6{0,l,2};i#A; 

where 7^^ = 7|\ + - e. 
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FIG. 9: Numerical fitting of |1) — >■ |2) spectral line (the same 
as traces shown in Fig. [S]) with different vales of e. The fol- 
lowing parameters (as in Table[T| are used for this simulation: 
Fio = 27r X 7MHz, = 2-k x llMHz, 7fo = 27r x 7MHz and 
7|'o = 27r X 16MHz. 



To summarize, the master equation including both 
energy relaxation and dephasing is (now back in the 
Schrodinger picture) 

P = -i[Hs, p]+ n (a'oipo-21 + O"i2pcrio) 
H — Y (^'^oi/oo'io - o-iip - pail ) 

"I Y (2'^12/°^21 - 0-22/3 - Pcr22) 

+ {2aiip(Jii - (Jiip - pail) 

, 720 



(2CT22PCr22 — Cr22P — pcr22) 
+ 2 (O'llPO-22 + a22P(Jll) , 



(B12) 



which we use to fit the experimental data with fixed de- 
phasing rates, relaxation rates and different values of e. 

As shown in Fig. 13 a finite value of e of about 27r x 
5MHz provides the best fitting (the (blue) dashed dotted 
curves). Values such as e = MHz and e = 27r x lOMHz 
provide reasonably good fittings, while larger values, such 
as £ = 27r X 15 MHz result in much worse fittings. We 
have also simulated the Autler-Townes spectra of Fig. [S] 
and we see no significant difference for values < e < 
2ti X lOMHz, while for larger values the fittings become 
worse. Changing the value of other parameters (such as 
"i2(y) does not produce better results either. We conclude 
that in our experiment £ is approximately 27r x 5MHz. 

Next, we suggest that the dephasing in this system can 
be understood as originating from two processes. Since 
the system is an oscillator with only weak anharmonic- 
ity, the first process is the fluctuation of the frequency 
LJQ of the oscillator , which depends on external fluctu- 
ating parameters such as the bias fiux $dc (as defined 
in Section |n|. The second source of dephasing are vir- 
tual transitions between the qubit and the environment 
(electromagnetic degrees of freedom, two-level systems, 
etc.). The first process can be experimentally controlled 
in a more straigthforward way, for example by filtering 
of the bias lines and by using magnetic-flux pinning sub- 
strates, while the second process is of a more fundamental 
nature. For both processes, we can use the general for- 
malism presented above to find the expressions for the 
dephasing Liouvilleans. 

The first process will be characterized by a single de- 
phasing parameter 7^91 since we note that 



^^10 (^) = l^^'J'dc(i) 

9(^10 +UJ21) 



dujo 



■(5$dc(i) « -56^21(0, 



SuJ2o{t) 



dc 



6^dcit) w 2Sujioit) w 2SL02iit), 
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where the derivatives are taken at the flux bias point 
and (5$cic(i) is the fluctuation in the dc-bias magnetic 
field (see also Section HT]). From Eqs. dH [BTl [B8l and 
E9l) we have = 2/1-2 jt dt' {dujo{t)6iOa{f)) and 







-7ho^oi -47ho^02 

-47hoP20 -7hoP2i 

(B13) 

Alternatively, we can start with the known expression for 
the dephasing of the harmonic oscillator— 

^Ho[/5] = ^ [2a^apa^a - {a^a)^p - p{a^af] , (B14) 

and truncate the creation and annihilation operators of 
the harmonic oscillator to the lowest three levels, 



(B15) 



In this three-level Hilbert space, Eq. (|B14I) is identical 










1 




1 


a = 


^2 




V2 








to Eq. (|B13|) . 

The second process can be characterized by two pa- 
rameters, and and, for processes involving virtual 



transitions with the environment the Liouvillean can be 
obtained again from the general expressions Eqs. (|B6[ 
IB7[ IB8[ and IB9p for the particular case of zero interlevel 
correlations. The resulting expression is well-known also 
from quantum optics^i^S, 



^virt[p] ^ J2 ^ (2c^iiPCrii - <^jjP - P<^jj) ■ 



(B16) 



With these assumptions in mind, we then have >Cdoph[p] 
Cko[p]+ ^virt[p], which results in jf^ 



47S0 

s = 47.'^ 



r^, and 7J1 



7S0 



^ 2 ■ 



7S0- 



rr,7ro 



This means that 



HO' 



and if we take e = 27r x 5 MHz we find 



2tt X 1.25 MHz, Tf = 
MHz. 



27r X 5.75 MHz, and 



-pip 
^ 2 



7ho 



27r X 11 



Appendix C: The steady-state density matrix 

To find the elements of the steady-state density ma- 
trix, it is possible to solve analytically (using Mathemat- 
ica) the equation p*^*'*) — for the case Ac = 0. The 
complete expressions are complicated, but they can be 
simplified and put in a form amenable to physical inter- 
pretation in certain limits. In this appendix we consider 
the approximation fl^ <C 7201 7ioi 7io720 < ^l- We find: 



J 



(st) 
PlO 



(st) 

P20 



(st) 

P21 



(st) 

Pll 



(st) 

P22 



(st) 

Poo 



(st)* 

Poi 



(st)* 

P02 



(st)* 

P12 



(2Ap 



-4A2 + 2zAp(7i 



720) 



1^? 



-4A2 + 2iAp(7io-f 72o) + ^^2' 
(4A2r2i(rio - 710) + 2iApTio^l 



r2i{Tio + 72oM)nl 



-4A2 



(4Ag7io + 720^^g) 
rio(-4A2 + f]2)2 ' 

(4A2(rio ^ 710) - (Eio -I- 72o)»g) nl 
rio(-4A2 + n2)2 

i-p[f-pif. 



(Cl) 
(C2) 

(C3) 

(C4) 

(C5) 
(C6) 



As a quick consistency check, one sees that the elements Inserting now the expressions of Eqs. (jCHIC6|) and using 



(st) 

the expressions Eq. (1551) and Eq. ([57]) if we restrict 
ourselves to 7201 7iO' 7io720 ^ ^c- Using the expressions 
above, we can calculate the fidelity Eq. (p9)) with respect 
to the dark state in the rotating frame. 



Eq. dSll) and p}^^' Eq. can be obtained from the approximations flp < 72O' 7iO' 7io720 < i* can 

be shown by algebraic calculations in Mathematica, that 
the fidelity is close to 1 for any detuning Ap. 

More insight into the structure of the state can be ob- 
tained, however, in the resonant case, Ap = 0. In this 
situation the density matrix elements Eqs. (|CHIC6|) be- 
come 



cos 29 ^ (st) 
o v^oo 



sin2e (st) 
o yP2o 



(st) 
Pll 



(C7) 



(st) 
PlO 



(st)* 

Poi 



»720^p 



(C8) 
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(st) 

P20 



(st) 

P21 



„(st) 

Pll 



(st) 

P22 

(st) 

Poo 



(st)* 

P02 



(st)* 

P12 = 



-ir2i(rio 



(rio + 72o)»g 



10 
(st) 

Pll 



(st) 

P22 • 



(C9) 

(CIO) 

(Cll) 

(C12) 
(C13) 



Now for an ideal dark state \D) in the limit ftp ^ ftc 
we have \D) = (1 - nl/2nl)\0) - (r2p/17c)|2) and the 
corresponding matrix elements are 







(C14) 




_ *{\D)) _ 


(C15) 






(C16) 



0, 
1 - 



(C17) 
(C18) 
(C19) 



A simple visual comparison between the two sets of ele- 
ments shows that they coincide up to first order in rip/flc 
(assuming 720/^^0 is of the same order or smaller). Even a 
better coincidence, up to second order in fip/fic, could in 
principle be reached provided that (72o/f^c) ^ i^p/^c)'^, 
(720/rio) < {np/n,), and (Fai/rJc) ~ (Tw/^c) ~ 
{Hp file). This agrees with the intuition that the effect 
introduces deviations from the ideal dark state second 
order in ilp/ilc are the decoherence effects associated 
with state |2) (dephasing and relaxation to state |1)). 
In conclusion, it is not surprising that states which are 
reasonable close to true dark states are produced in this 
experiment. 
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